One-dimensional conductance through an arbitrary potential 
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The finite-size Tomonaga-Luttinger Hamiltonian with an arbitrary potential is mapped onto a 
non-interacting Fermi gas with renormalized potential. This is done by means of flow equations 
for Hamiltonians and is valid for small electron-electron interaction. This method also yields an 
alternative bosonization formula for the transformed field operator which makes no use of Klein 
factors. The two-terminal conductance can then be evaluated using the Landauer formula. We 
obtain similar results for infinite systems at finite temperature by identifying the flow parameter 
with the inverse squared temperature in the asymptotic regime. 
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The conductance of one-dimensional (ID) conduc- 
tors at low temperatures is determined by a sophisti- 
cated interplay between impurity scattering and electron- 
electron interaction. The question was first raised about 
thirty years ago by Luther and Pcschei^ and Mattis 2 
in the context of enhanced conductivity, present in ID 
structures at low temperatures. Their results were later 
confirmed and extended using renormalization group 
techniques i&4 Recently, evidence was found that the 
conductance in single-walled carbon nanotubes 5 - and 
GaAs/AlGaAs quantum wires 6 - behaves according to the 
so-called Luttinger liquid theory^ 

The above theoretical descriptions were only valid for 
large systems at low temperatures due to the perturba- 
tive treatment in the impurity strength. This shortcom- 
ing was overcome by Monte-Carlo calculations^ and ap- 
plying the thermodynamic Bethe ansatziS Nevertheless, 
the latter method already fails for a double impurity. 

To describe general external potentials, it thus seems 
favorable to tackle the problem from the other end, i.e., 
to treat the impurity exactly but include the electron- 
electron interaction perturbatively. Since the interaction 
of ID systems is marginal in the renormalization group 
sense, 10 this has to be clone with care. One approach is 
to use scattering theory within the Born approximation 
and renormalization group techniques The interac- 
tion can also be included via exact renormalization group 
equations based on the one-particle irreducible effective 
action H 

Recently, we have employed a flow equation scheme 
for Hamiltonians 14 which is applicable for small electron- 
electron interaction but arbitrary impurity strength. 15 
For infinite systems at zero temperature, the approach 
yielded a logarithmic expansion for the spectral weight 
at the impurity site involving the exact boundary expo- 
nent. In this Brief Report, we will extend the analysis 
and recover the exact power-law behavior. For this, we 
will derive the exact Green's function of the Tomonaga- 
Luttinger model without making use of the bosonization 
formula. We will also demonstrate how the conductance 
can be obtained as a function of the system size and the 
temperature, respectively, using the Landauer formula. 16 

We start our discussion from the spinless Tomonaga- 



Luttinger Hamiltonian of finite length L 
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where vf denotes the Fermi velocity and v q — V— q is 
the Fourier transform of the electron-electron interaction. 
The sum is over all non-zero wave numbers q ^ and the 
normalized operators of the electron density fluctuations 
b q ^ obey canonical commutation relations. They are 
given by 6 ± | g | = -±= J2 k 4,± c fct|g|,± where c% ± are the 
fcrmionic annihilation (creation) operators of the plane 
wave with wave number k of the right- and left-movers 
and k = 2im/ L is unbounded (n € Z). 

We now perform a continuous unitary transformation, 
U(£), such that the interaction Hamiltonian H ee tends 
to zero for £ — > oo. In differential form, this transforma- 
tion can be characterized by the so-called flow equations 
dgH(£) — [r)(l), H(£)] with the anti-Hermitian generator 
7](£) = — ry(£)i^ The parameters of the Hamiltonian |QJ 
thus become functions of the flow parameter £. This will 
be indicated by dropping the upper zero which was in- 
troduced in order to specify the initial conditions. 

In this work, we choose the generator as rj = [Ho, H ee ] 
even in the presence of an external potential^ This 
yields 



(2) 



with i] q (£) = -2u q (£)uq(£) = -sgn(u g )^[siiili(4^£ + 
Cq)]- 1 where smh(C ? ) = u 2 q /2uj G q \u a q \ and Q q = u g (£ = 
oo). The fixed point dispersion uj q follows from the flow 
invariant uj q (£) — u 2 q {£) — v 2 F q 2 {l + v q /irvF) as u q {£) — > 
for £ — * oo. For attractive electron-electron interaction 
we impose the lower bound v q > —ttvf in order to avoid 
phase separation^ 
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The ^-dependent parameters then read 

LJ q (£) = c q (£)oj° q + Sq (£)u° q , u q {£) = c q (£)u° q + s q (£)cu° q 

(3) 

where c q (£) = cosh(E q (£)) and s q {£) = s'mh(E q (£)) with 

/ tanh(2w^)if2 + i 



E q {i) = / dHr} q {£!) = In 



Vtanh(2d^)x- 



(4) 



and the Luttinger liquid parameter K q = (uj q — u q )/(uj q + 

In order to evaluate correlation functions, the observ- 
able has to be transformed by the same continuous uni- 
tary transformation, i.e., diip±{x,£) = [??(£), ip± i x , £)]• 
Starting from the fcrmionic representation of the field 
operator, i.e., ip±(x) = L -1 / 2 J2k elkx °k,±, the flow equa- 
tions generate an infinite hierarchy of operators. Never- 
theless, this expansion factorizes and only involves one 
g-dependent function of £, ip^ (£) . The expansion can be 
summed up to an exponential and is given by 



ip±(x,£) = ip±(x)exp [ 



(•5) 



where we defined <j> q {x) = e zqx b q — e %qx b q . Notice that 
[tp±{x), 4> q (x)] = 0. The flow equation for ^{£) reads 



d i( p q L (£) = - Vq (£)(Q(Tq) + tp± q (l)) 



(6) 



where Q(x) denotes the Heaviside step function. One 
obtains the following solution 



<pf(e)=e(±q)(c q (£)-l)-e(Tq)s q (£) 



(7) 



with c q {£) and s q {£) defined below Eq. ©. Notice that 
no Klein factor is needed in the above representation, Eq. 
©■ 

We now define normal ordering : ... : such that 
the bosonic operators b q (b q ) stand left (right) from the 

fermionic operators tp± (x). With the Baker- Hausdorff 
formula e A+B = e A e B e -\ A > B \l 2 anc j the relation 

*±(x)f(b\) = f(b\ + Q{±q)n q l l 2 e^)i, ± {x) , 
which holds for an arbitrary function /, we obtain 



(8) 



ip±(x,£) = :ip±(x,£): exp [- 



q>0 



(9) 



With this normal ordering procedure, we can thus ex- 
tract the anomalous scaling factor L~ s W, which follows 
in the limit L — ► oo from a constant potential with fi- 
nite cutoff g c , i.e., s{£) = s q (£) for q < q c and zero oth- 
erwise. From the scaling behavior of the transformed 
field ip±{x,£ = oo) - known from bosonization - we infer 
that the normal ordered field operator scales like a free 



fermionic field. This will enable us to substitute the nor- 
mal ordered transformed field operator by the bare field 
for weak electron-electron interactions. 

To evaluate the Green's function, we need to normal 
order the operator 

g£ 0, x';t) = i/>± (x', 0,£ = oo)V>± (x, t,£=oo), (10) 

where the time dependence is defined in the Heisen- 
berg picture with the fixed point Hamiltonian H* = 
J2 q & q b q b q . We obtain 

g£(x,x';t) = :g£(x,x';t): g(x-x',t) , (11) 
with (s q = s q {£ = oo)) 

g(x, t) = exp [ - J2 T^C 1 " cos^arje* 5 '*)] . (12) 



<j>o q 



The Green's function G±(x,x' ,t) 
thus reads 



-i{g£{x,x';t)) 



G^(x,x' ) t) = -i(^ ± {x',0)ij ± (x,t))g(x-x',t) (13) 



where 
E Uqbib q 



denotes the ground-state average of H* = 



■j q - J q w q u q- 

Up to now, we have assumed a finite interaction cut- 
off q c which yields well-defined expressions. For con- 
stant interaction potential v q — vQ(q c — \q\) and q c — > 
oo , the remaining expectation value can be evaluated 
using the Kronig relation which relates H* to a free 
fcrmionic Hamiltonian with linear dispersion. 20 With 
u q = v c \q\, v c denoting the charge velocity, and Ck,±{t) = 
Ck,±e T1Vc( - kTkF ^ , we obtain for L — » oo 

„±ikp (x— x') 



(4 (x',0)i> ± (x,t)) = 



—i 



2n ±(x — x') — v c t — iO 



(14) 



This yields the well-known result for the Green's func- 
tion of the Tomonaga-Luttinger model. 

To complete the discussion on the observable flow, 
we also transform the ladder operators Ck,± according 
to diCk,±{£) — [rj(£),Ck,±(£]]- Again, an infinite hierar- 
chy of operators is generated which is only character- 
ized by one g-dependent function of £, ip^ {£) . Defining 
4>k', q = b q S k ',k~q - b q d k >,k+q, the series is given by 



oo 1 n 



(15) 



with fco = k. The q-dependent function of £ is again given 
by Eq. 0. Notice that $±(1) = L~ 1/2 £ fe e ifc:E c fe , ± (f). 
For a step function potential v q = vQ(q c — q) in the limit 
L — > oo, we therefore have 



c kt± (£) = :c kt± (£): 
where s(£) as defined below Eq. JJjJ. 



(16) 
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We will now add an arbitrary external potential to the 
initial Hamiltonian, Eq. ft|. The additional contribu- 
tion H e shall consist of a forward and a backward scat- 
tering part, i.e., H e = Hp + Hb- The forward scatter- 
ing contribution can be expressed by the bosonic density 
fluctuations )Stt 

Hr = ? nV=L-/^» W («<-^ + e - 6i ). (17) 



It is interesting to analyze the asymptotic behavior of 
Ub{x, £) for L — ► oo and I — ► oo. For this, we will choose 
a step function potential v q — vQ(q c — q) with finite 
interaction cutoff q c and also drop the spatial component 
in the following. Considering first the limit L — ► oo, 
the finite interaction cutoff q c — 27m c /L shall remain 
constant, such that n c is proportional to the system size 
L. With X)"=i V n — >Tnn c + C for n c — > oo, where C is 
Euler's constant, we readily find (X = iC g for < q c ) 



The backward scattering contribution reads 



Hp 



dxU B (x)(ipUx)ip_(x) + h.c), (18) 



with the right- and left-moving field, ip±(x). 

The flow equations for Hp close and the differential 
equations for the forward scattering potential read 



d i U?(x,t) = r h (£)UEJx,t) 



(19) 



We thus have U^{x,£ = oo) = ^/K q U Ffl {x) such that 
the forward scattering amplitude is enhanced for attrac- 
tive and reduced for repulsive electron-electron interac- 
tion. 

The flow of the backward scattering operator 



o ±2kF (xj) = iPi(x,iw T (x,e) 



(20) 



is determined by the flow equation of the field operator 
ip±(x,£), outlined above. Defining the ^-dependent Lut- 
tinger liquid parameter K q {£) = exp(2E q (£)), we have 



0±2k F (x > £) = ±2 k F (x) 



(21) 



L q V n i 



with cf> q (x) = e lqx b q — e lqx b q . Normal ordering the 
above expression such that b q (b q ) stands left (right) from 
O ±2kF {x ) £ = 0) yields 



±2kF (x 1 £) = :0 ±2kF {x,£): 

1-K q (£) 



(22) 



x exp 



q>0 



Above, it is shown that the normal ordered trans- 
formed field operator scales like a free fermion field. 
For small electron-electron interaction ^jK q {£) — 1 — > 0, 
the normal ordered backscattering operator can thus be 
approximated by the bare backscattering operator, i.e., 
:0 ±2kF (x,£):^0 ±2kF (x). 

This approximation provides a suitable truncation 
scheme for the Hamiltonian flow. We can then define 
an effective i'-dependent backscattering potential as 



U B (x,£) = C/ S (x)exp [£- 

q>0 



K q {£) 



Ub oc L 



l-K 



(24) 



Notice that the renormalization of the external potential 
depends linearly on the interaction strength v. 

For an infinite system in the limit £ —> oo, the sum in 
the exponent of Eq. (|23|l turns into an integral accord- 
ing to J2 q — > jjp J dq. For a constant electron-electron 
interaction with sharp cutoff q Cl we then obtain 



U B oc £^- k )' 2 



(25) 



Notice that the asymptotic limit £ — > oo and the ther- 
modynamic limit L — > oo interchange. The commutation 
of the two limits does usually not hold even if one consid- 
ers impurity properties. But in contrary to the treatment 
of, e.g., the spin-boson model of Ref. |2J, here we do not 
neglect finite-size corrections in the flow equations. 

For repulsive electron-electron interaction, an initial 
delta-impurity strength thus scales to infinity for L — > oo 
and T — » 0. The ground-state of the fixed point 
Hamiltonian is then given by the ground-state of an 
electron gas with fixed or open boundary condition. 15 
Concerning the operator flow, this change in boundary 
conditions can be accounted for by the substitution 



b-\ q \ — > b\ q \ in the operator (f> q (x), defined below Eq. 
(JSJ), i.e., there are only right-movers. With the same 
normal ordering procedure as outlined for the case of 
periodic boundary conditions, this yields the well-known 
Green's function of the Tomonaga-Luttinger model with 
open boundaries^ We thus obtain the exact power law 
behavior for the spectral weight at the impurity site. 

So far, we have mapped the TL model with external 
potential onto a non-interacting bosonic bath with dis- 
persion relation ui q = vpq^Jl + v q /irvF and a renormal- 
ized external potential. The mapping holds in the limit 
of weak electron-electron interaction. Linearizing the dis- 
persion relation around q — 0, one can map the non- 
interacting bosonic bath onto a non-interacting Fermi 
gasmQ Notice that, up to terms linear in v q and approxi- 
mating :c k ,±{£) ■ ~ Ck,±i these are the same fermions we 
started with, see Eq. 0fy M 

The conductance G is now given by the Landauer for- 
mula for a one-channel, spinless electron system at T = 0, 

2 

G = \T i where T denotes the transmission probabil- 
ity through the renormalized external potential>iSi24 By 

(23) 

this, we assume that the coupling of the interacting sys- 
tem to the non-interacting leads does not contribute to 
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the transmission probability. T is thus entirely deter- 
mined by the renormalized potential which corresponds 
to perfect coupling. 25 Nonadiabatical coupling can be dis- 
cussed by considering different boundary conditions of 
the fermionic field operator. 

One can also obtain results at finite temperatures if 
one identifies the asymptotic flow parameter £ with the 
inverse temperature squared, i.e., 

Vlocl/T . (26) 

For this, it is crucial to choose the generator r\ "canon- 
ically", i.e., based on the commutator of a diagonal 
and a non-diagonal part. Then, in the case of a sim- 
ple model, the flow parameter £ is related to the in- 
verse squared energy uncertainty of the basis states of 
the fixed point Hamiltonian^4 This interpretation is con- 
firmed in case of the spin-boson model as the effective 
tunnel-matrix element A* is asymptotically reached ac- 
cording to A(£) — A* oc 1/V«i2i Furthermore, it turns out 
that spectral functions, initially depending on the energy 
lo and on the flow parameter £, can be characterized by 
functions of the scale-independent variable y — toVJ for 
£ — ► oo. 26 This also holds for the Tomonaga-Luttinger 
model, see Eq. (gj. 

Physical quantities should not depend on energy shifts 
of the spectrum which are smaller than the energy scale 
given by the finite temperature. Therefore, the flow equa- 
tions do not have to be integrated up to £ = oo, but 
it suffices to halt the integration at a finite £. Since 
the electron-electron interaction is exponentially small 
for large £ and not marginally relevant at finite tempera- 
tures, we can neglect it. This leaves us with a one-particle 



problem with ^-dependent parameters. The temperature 
is introduced by substituting £ by T~ 2 . 

The above procedure yields the same result for the 
conductance as obtained by Ref. in case of a delta 
impurity, but is also applicable for general potentials. 
From Eqs. JUJ, (JS5J, and (JJSJ, we can also infer the 
finite-size scaling relation L oc 1/T. 

Finally, we want to mention that Refs. ^3 and 111 ob- 
tain an effective non-interacting model with an oscillating 
and bounded effective potential. But whereas these ap- 
proaches depart from a non-interacting system with an 
impurity and then include the electron-electron interac- 
tion, we start from an interacting system and add the 
impurity later, which might explain the different fixed- 
point Hamiltonians^i 

To conclude, we have used Wegner's flow equation 
method to map the Tomonaga-Luttinger Hamiltonian 
with an external potential onto a non-interacting Fermi 
gas with an renormalized external potential. The map- 
ping is valid for weak electron-electron interaction and 
is based on an alternative bosonization formula for the 
transformed field operator derived by the flow equa- 
tion method. Calculating the transmission amplitude 
from the non-interacting representation, the Landaucr 
formula then provides an analytic expression for the one- 
dimensional two-terminal conductance for finite size sys- 
tems. Interpreting the flow parameter in the asymptotic 
regime as inverse temperature squared yieldes an analytic 
expression of the conductance also at finite temperatures. 
The formalism applies to arbitrary external potential. 

The author wishes to thank V. Meden, A. Mielke, and 
F. Wegner for helpful discussions. 
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